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Non Gaussian states and processes are useful resources in quantum information with continuous 
variables. An experimentally accessible criterion has been proposed to measure the degree of non 
Gaussianity of quantum states, based on the conditional entropy of the state with a Gaussian 
reference. Here we adopt such criterion to characterise an important class of non classical states, 
single-photon added coherent states. Our studies demonstrate the reliability and sensitivity of this 
measure, and use it to quantify how detrimental is the role of experimental imperfections in our 
realisation. 



Quantum information offers a different viewpoint on 
fundamental aspects of quantum mechanics: it aims to 
assess and exploit the quantum properties of a physical 
system as a resource for different, and wishfully more 
efficient, treatment of information. Indeed, within the 
framework of quantum information with continuous vari- 
ables pQ, nonclassical states of the radiation field repre- 
sent a resource and much attention has been devoted to 
their generation schemes, which usually involves nonlin- 
ear interaction in optically active media. 

On the other hand, the reduction postulate provides 
an alternative mechanism to achieve effective nonlinear 
dynamics; if a measurement is performed on a portion 
of a composite entangled system, the other component 
is conditionally reduced according to the outcome of the 
measurement. The resulting dynamics may be highly 
nonlinear, and may produce quantum states that can- 
not be generated by currently achievable nonlinear pro- 
cesses. Conditional measurements have been exploited to 
engineer nonclassical states and, in particular, have been 
recently employed to obtain non-Gaussian states. 

While Gaussian states, defined as those states with a 
Gaussian Wigner function, are known to provide useful 
resources for tasks such as teleportation [2j [3], cloning 
[1H6], or dense coding [7H2], there is an ongoing effort 
to study which protocols are allowed by non- Gaussian 
resources. The most notable example is certainly their 
use for an optical quantum computer [TO] [11] , alongside 
with their employment for improving teleportation [T2l - 
[14], cloning [15] , and storage [16]. Several realisations of 
non- Gaussian states have been reported so far, in partic- 
ular from squeezed light [T7H25] , close-to-threshold para- 
metric oscillators [26] [27] , in optical cavities [28] , and in 
superconducting circuits [29]. Non-Gaussian operations 
are also interesting for tasks as entanglement distillation 
[30] [31], and noiseless amplification [32j[33] which also are 
obtained in a conditional fashion, accepting only those 
events heralded by a measurement result. 

In principle, non- Gaussianity is not directly related 



to the non-classical character of a quantum state and, 
in turn, classical non-Gaussian state may be prepared, 
e.g. by phase-diffusion of coherent states or photon sub- 
traction on thermal states [34]. On the other hand, 
in the applications mentioned above it is the presence 
of both non-Gaussianity and non-Classicality which al- 
lows for enhancement of performances. Therefore, de- 
Gaussification protocols of interest for quantum informa- 
tion are those providing non-Gaussianity in conjunction 
with nonClassicality. 

In this work we address the conditional dynamics in- 




FIG. 1: Layout of the experiment. An optical parametric am- 
plifier (OPA) is injected with a coherent state of variable am- 
plitude || a || in the range [0, 1.5]. This is realised by a 100/xra 
thick slab of potassium niobate, pumped by a frequency- 
doubled mode-locked TkSapphire laser (A p =425nm, pulse du- 
ration 230fs, repetition rate 800kHz). Our OPA is driven in 
frequency-degenerate and non-collinear regime, so to gener- 
ate an idler at the same wavelength A = 2A P as the coherent 
seed; this is then spatially filtered with a single-mode fibre, 
spectrally filtered by a diffraction grating and a slit, indi- 
cated as F in the figure. Finally, the idler is detected by an 
avalanche photodiode APD. The observation of the output 
conditioned by an APD count results in single-photon addi- 
tion. The quantum state of the output is reconstructed by 
homodyne detection (HD). Mode-matching with the local os- 
cillator (LO) employs polarisation: the signal and the LO are 
first matched on a polarisation beam splitter, and then com- 
bined using a half-wave plate and a second polariser so to 
realise an accurate 50:50 intensity splitting. 
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duced by the so-called photon addition as a protocol 
to generate nonclassical non-Gaussian states. We quan- 
tify experimentally the amount of non-Gaussianity ob- 
tained by adding a photon to a coherent state [191 [351 
[37] . Differently from previous investigations [35j l38l 
[41], we can explicitly address the two aspects of non- 
Gaussianity and non-Classicality at once. For the for- 
mer, we adopt the non-Gaussianity measure S[g] pro- 
posed in [42] [43], defined as the quantum relative en- 
tropy between the quantum state itself g and a reference 
Gaussian state r having the same covariance matrix as g. 
Given this choice of the reference Gaussian state, we have 
that Tr[^>logr] = Tr[rlogr], as logr is a polynomial of 
order at most two in the canonical variables [42j [44] . We 
thus find 



5[ e ]=S(Q\\T) = Tr[ Q QagQ-]QgT)] 
= S(t) - S(q) 



(1) 



that is, 5[q] is simply equal to the difference between 
the von Neumann entropy of r and the von Neumann 
entropy of g. In Ref. [42] it has been shown that this 
measure is non zero only for non-Gaussian states. It is 
also additive under tensor product, invariant under uni- 
tary Gaussian operations, and in general it does not in- 
crease under generic completely positive Gaussian chan- 
nels. This measure is somehow preferable to that based 
on the Hilbert- Schmidt distance [45 in a quantum in- 
formation context, since it is based on an information- 
related quantity. We note, however, that a mixture (e.g. 
doubly peaked) of classical states can also be strongly 
non-Gaussian: we therefore adopt an additional "non- 
classicality " criterion. 

Several measures of non-classicality have been pro- 
posed in literature [46 -49 , for our purposes we consider 
as a witness a quantity v[g] related to the negativity of 
the Wigner function. This is normalised to a reference, 
which we choose to be a single photon state W\{x,p). 
The non-classicality is then defined as 



min (W(x,p)) 
min (Wi(x,p)) 



(2) 



This reference has been chosen since it has the lowest 
value within the class of states we consider. While this 
does not constitute a measure, it acts as a witness for 
non-Classical states whenever u[g] > 0. The choice of us- 
ing the single-photon as a reference is dictated by the 
need of a measure which does not depend on the con- 
vention for the quadratures. Moreover, it sets to unity 
the highest value of v[g] attainable in the class of states 
under investigation. 

A conceptual scheme of the experiment is shown in Fig. 
1: an input coherent beam \a) is injected in an optical 
parametric amplifier (OPA). This is a three- wave non- 
linear interaction between a pump beam and the input 
beam (usually called the signal s) which results in the 
generation of a third beam called idler (i). When the 
pump is an intense beam, we can treat it as a classical 



field: the output state of the s and i modes can then be 
expressed as the application of the squeezing operator: 



£ s ,i(r)=exp (r(a\a\ - a^a s )) 



(3) 



to the input |a) s |0)$. Here, r is the squeezing param- 
eter, which depends on the pump intensity, the crys- 
tal length and its non-linear coefficients, and attained 
a value r~0.105 in our experiment; we can then approx- 
imate S Si i(r) taking the limit of weak nonlinearity: 



S Syi (r)~I + r(ala\) - r(a i a s ). 



(4) 



We now put an avalanche photodiode (APD) on the 
idler beam, and accept only those events when a click is 
registered. Since the idler was originally in the vacuum 
state, the only term which can give a contribution in 
S Si i(r) is the second one. Therefore, the detection of a 
single photon on the idler heralds the addition of a single 
photon to the coherent state, transforming it into 
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(5) 
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FIG. 2: Non-Gaussianity S[g] (upper panel) and non- 
classicality is[g] (lower panel) as a function of the ampli- 
tude \a\ of the input coherent state for different values of 
the squeezing parameter r (dashed lines); from top to bot- 
tom r = {0.15,0.30,0.45}. The black solid line corresponds 
to the non-Gaussianity of the ideal photon added coherent 
state, that is to the limit r — »• 0. 
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FIG. 3: Experimental Wigner functions for increasing values of a. The output states are reconstructed by a maximal likelihood 
algorithm [51] which interpolates 800,000 data points sorted according to their phase into 12 histograms. The effectiveness of the 
sorting algorithm sets a lower bound \a\ ~ 0.5, so that the oscillations due to interference are much larger than low-frequency 
noise fluctuations. Notice that for ||a||=0 this noise can be compensated by using a moving average technique. 



in the ideal case. In practice, we need a careful analysis of 
those processes which spoil the photon addition and the 
non-Gaussianity of the resulting state. Here, we follow 
closely the model presented in Refs.[20 | 123 } 130], 

The detection on the idler beam is performed by an 
APD that can not resolve photon number. In the limit 
of small detection efficiency, we can approximate the de- 
tection process as the application of the ai annihilation 
operator on the idler mode. This is actually the case in 
our experiment, where the overall detection efficiency is 
less than 10%, due to spatial filtering (<75%), spectral 
filtering (<30%), and limited efficiency of the photodiode 
(55%). 

In Fig. [2] we plot S[g] and v[q] as a function of the 
coherent amplitude a, for different values of r. We ob- 
serve that the two trends resemble closely, suggesting 
that the non-Gaussianity induced by photon-addition is 
essentially non-classical and thus useful for quantum in- 
formation processing. It can be also observed how both 
non-Gaussianity and non-classicality decrease by increas- 
ing the squeezing parameter; this can be explained by ob- 
serving that, as shown in Eq. Q, for low values of r, the 
squeezing operator adds only one photon on each arm, 
while by increasing r we have to consider also the possi- 
ble addition of many photons. Due to the lack of photon 
number resolution, the detection will be affected by the 
presence of higher-number emission from the squeezing 
process, eq. In this case, conditioned on a click from 
the idler beam, the signal will be in a highly mixed and 
thus less non- Gaussian and also non-classical state. In 
the ideal limit of r — » the non-Gaussianity of the state 
is exactly equal to the one of the ideal photon added 
coherent state in Eq. ([5]). However, this goes to the ex- 
penses of the success rate, and a compromise between 
non-Gaussianity and count rate has to be found. 

Beside the role played by the squeezing, we have to take 
into account the other imperfections that are present in 



our experimental setup. In the OPA there might occur 
a certain modal mismatch between the pump and the 
input: this results in a parasitic amplification that intro- 
duces excess noise on the signal and idler modes. The 
process is modelled as a non-degenerate OPA driven at 
a weaker strength 77*, where < 7 < 1 and 7 ~ 0.425 
in our experiment. The amplification couples the modes 
s and i with two other modes s' and z', initially in the 
vacuum state. The complete description takes the form: 



S s , i (r)5 s , i ,( 7 r)5 sV ( 7 r)|a) s |0) i |0) s H0) i 



(6) 



The parasite modes s f , i' are not observed in the experi- 
ment, therefore we have to trace over them to obtain the 
output density matrix. 

Accurate spatial and spectral filtering is performed so 
that the mode detected by the APD is matched with 
the input mode, that we detect by balanced homodyne; 
nevertheless, this task can be accomplished only with a 
limited efficiency £ which in our setup takes the value 
£ ~ 0.96. In formulae, we will have an output state g s y 
on the signal mode when the trigger count came from 
the correct mode, and a different state g SjX heralded by 
a faulty trigger event; the overall state is: 



(7) 



Notice that dark count rates from the APD play a negli- 
gible role (~ 10 counts/s with an overall rate ~ 1—4- 10 3 
counts/s), as we used a gated detection triggered by the 
cavity dumping electronics of our laser. Homodyne de- 
tection has a limited efficiency as well, coming from op- 
tical loss, non-unit detector yield, and mode-matching 
between the local oscillator and the signal. The overall 
efficiency is, in our case, of the order of 77 ~ 0.71; this is 
modelled as transmission through a beam splitter with 
transmittivity t 2 =rj. Examples of the measured Wigner 
quasi-distributions are illustrated in Fig. [3] 
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In Fig. [4] we plot S[g] at fixed values of the coherent 
state amplitude a = 0.5 and of the squeezing param- 
eter r = 0.15 as a function, respectively, of the noise 
parameters 7, £ and 77 chosen in ranges relevant for our 
experimental setup. We observe as expected that 5[g] 
decreases monotonically with 7, while it increases mono- 
tonically with £ and 77. For the values that characterize 
our experiment, the homodyne efficiency 77 is the source 
of imperfection that affects in the most detrimental way 
the non-Gaussianity of our states. 

Finally we evaluated S[g] from the experimentally re- 
constructed output states for different coherent state am- 
plitudes \a\; the results are shown in the upper panel of 
Fig. [5] . The dashed line shows the description provided 
by our model when taking into account all the noise pro- 
cesses described above. The values of parameters used for 
the curve are obtained from a fit of the experimental data: 
r = 0.105, 7 = 0.425, £ = 0.96 and 77 = 0.71. The average 
fidelity between the reconstructed and the modelled state 
is 0.989±0.006 [50]. Concerning the non-Gaussianity, the 
agreement between the experimental data and the model 
is satisfactory, and we can observe, as expected, a de- 
crease as the input intensity \a\ increases. As shown, the 
effect of the single-photon addition is more relevant for 
quantum states with a small number of photons, and be- 
comes only a small perturbation for higher average pho- 
ton number. 

In the lower panel of Fig. [5] we observe the behaviour 
of v[q] as a function of the amplitude a. The exper- 
imental results confirm that the two quantities, non- 
Gaussianity and non-classicality have a similar behaviour 
and then that the non-Gaussianity induced by this 
photon-addition operation is essentially non-classical. As 
a general remark, we notice that the logarithmic term in 
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FIG. 4: Non-Gaussianity S[g] as a function of the noise 
parameters of the experimental setup for fixed amplitude 
I a I = 0.5 and squeezing parameter r = 0.15. In details: the 
blue dot-dashed line corresponds to S[g] as a function of 7, the 
green dotted line as a function of £ and the red dashed line as 
a function of 77. The solid lines refer to the non-Gaussianity of 
the ideal photon-added coherent state with \a\ = 0.5 (upper 
black line), and to the non-Gaussianity of the state obtained 
by considering \a\ = 0.5, squeezing parameter r = 0.15 and 
no imperfections (lower grey line). 
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FIG. 5: (Top): Non-Gaussianity S[g] as a function of the 
amplitude \a\ of the input coherent state. (Bottom): non- 
classicality is[g] - related to the minimum value of the Wigner 
function of g - as a function of a. The red points are the ex- 
perimental values from the reconstructed matrices. The black 
dashed line is obtained from our model including the main ex- 
perimental imperfections of our realisation. The parameters 
are chosen in such a way to fit the data of the non-Gaussianity: 
r 0.105, 7 0.425, f 0.96 and 77 0.71. 



the expression ([!]) amplifies the effect of small discrepan- 
cies with the model we present. This qualifies our mea- 
sure as a very sensitive one in those contexts where a 
good estimation of information resources is needed. In 
any case, our model is able to capture the essential fea- 
tures of the process, and provides us a tool to quantify 
how detrimental the imperfections are for the generation 
of these non-Gaussian resources. 

In summary, these experiments on single-photon added 
coherent states demonstrate the relevance of this recently 
proposed measure of the non-Gaussianity. This mea- 
sure appears as a reliable and sensitive way to quantify- 
ing experimental imperfections of de-Gaussification ex- 
periments. It furthermore allows to exhibit a link be- 
tween non-Gaussianity and non-Classicality in such ex- 
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periments. More generally, it would be useful to have 
at disposal a quantity able to capture both features, 
non-Gaussianity together with non-classicality, for any 
generic quantum states. Work along these lines is in 
progress and results will be developed elsewhere. 
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